Renormalization group study of one-dimensional systems with roughening transitions 
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A recently introduced real space renormalization group technique, developed for the analysis of 
processes in the Kardar-Parisi-Zhang universality class, is generalized and tested by applying it to 
a different family of surface growth processes. In particular, we consider a growth model exhibiting 
a rich phenomenology even in one dimension. It has four different phases and a directed percolation 
related roughening transition. The renormalization method reproduces extremely well all the phase 
diagram, the roughness exponents in all the phases and the separatrix among them. This proves 
the versatility of the method and elucidates interesting physical mechanisms. 
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I. INTRODUCTION 

Surfaces and interfaces may grow in a smooth way or 
alternatively in a rough fashion. The study of the physi- 
cal mechanisms originating these different behaviors has 
been the focus of an overwhelming number of recent stud- 
ies [0-|. The Kardar-Parisi-Zhang (KPZ) [|| equation is 
the minimal continuous model capturing the aforemen- 
tioned physics. In dimensions larger than two it may ex- 
hibit two different phases: a flat and a rough one. Sep- 
arating both of them there is a roughening transition. 
Apart of being a milestone in surface growth theory, the 
KPZ equation is also related to other interesting phys- 
ical problems: The Burgers equation in turbulence pi, 
directed polymers in random media [p| , and systems with 
multiplicative noise [Q among others. 

While the physics of the KPZ flat phase is very well 
understood, elucidating the properties of the rough phase 
has proven a puzzling problem ll|,|g,|9|]. In fact, standard 
field theoretical analysis finds an unavoidable difficulty: 
The rough phase regime is controlled by a strong cou- 
pling fixed point, not accessible by standard perturbation 
techniques. Therefore, from a field theory point of view 
not much can be concluded about the rough phase (re- 
cent interesting results in this direction can be found in 
P-[ll|). An alternative strategy has been recently pro- 
posed to deal with this elusive problem; namely a real 
space renormalization group (RSRG) approach. Its non- 
perturbative nature permits a direct access to the strong 
coupling regime and, in particular gives estimations of 
the roughness exponent in dimensions ranging from d = \ 
to d = 9 [O-O (see also flq] with results in very good 
agreement with the best numerical measurements [16|. 
Moreover, the same method has provided analytical evi- 
dence for the absence of an upper critical dimension for 
the KPZ strong coupling phase ||l3], which has been a 
highly debated subject |1^]. Additionally, it has also 



permitted to analyze the behavior of the simpler linear 
surface growth model, i.e. the Edwards- Wilkinson (EW) 
equation |14|. 

In this paper we intend to go further in the applica- 
tion and understanding of this new RSRG approach. In 
particular, we study a class of systems exhibiting rough- 
ening transitions even in one dimension, by renormalizing 
them with the new RSRG approach. Our motivation for 
that is twofold; on one hand we want to test the RSRG 
method (which was specifically devised to deal with KPZ 
growth) when generalized and applied to other physical 
situations, i.e. we intend to analyze its versatility to deal 
with different physical problems. On the other hand, by 
doing so we will perform a renormalization of the class 
of systems exhibiting a roughening transition in d = 1, 
that allows us to get some insight into their interesting 
physics. 

The paper is structured as follows. In section I we 
present the family of models exhibiting a one-dimensional 
roughening transition, and review their main properties. 
In section II we present a generic two-parametric model 
in this class suitable to be renormalized using the RSRG 
approach, and discuss in detail all the different phases 
and physical behaviors. In section III we briefly present 
the main traits of the RSRG approach, discuss its ap- 
plication to our model and present a detailed discussion 
of the results. In section IV the conclusions are pre- 
sented. Finally, in the appendix we study the connection 
between the one-dimensional roughening transitions dis- 
cussed previously and directed percolation. 



II. SYSTEMS WITH ONE-DIMENSIONAL 
ROUGHENING TRANSITIONS. 

In this section we review the class of models exhibit- 
ing a one-dimensional roughening transition. The firstly 
studied model in this class is the so called Polynuclear 
growth model analyzed by Kerstez and Wolf (KW) p8| . 
Its dynamics is defined by two successive steps: In the 



first one, particles are (parallelly) deposited with proba- 
bility p at each site of a one-dimensional lattice. In the 
second one, the terraces (kinks) formed by the previous 
deposition process, grow laterally in a deterministic way 
by u units (or less if less space is available on the ter- 
race). This is, kinks move deterministically increasing 
always the averaged height ||l^ . These two processes are 
iterated in time. For every u there is a critical value of 
p, Pc, such that for p > pc the surface grows homoge- 
neously, it is flat, and moves with maximal velocity p8[ . 
The roughness exponent in this flat regime is a = 0, and 
(3 = ||l^. On the other hand, for p < pc steps are 
less likely to annihilate and the surface becomes rough. 
In this phase the roughness exponent is estimated to be 
a ~ 0.5, and (3 « 0.33 compatible with their correspond- 
ing KPZ values. The roughening transition, as we discuss 
below, and as first pointed out by KW, is related to di- 
rected percolation (DP) pO|| . 

After the seminal work by KW other models proposed 
for rather different physical problems have appeared in 
the literature, exhibiting similar phenomenology. 

Alon et al. |2^] proposed a model with absorption of 
particles and desorption at the edges of grown islands. 
This mimics the fact that, in crystal growth, particles ab- 
sorbed in the interior of grown islands are more strongly 
bounded than particles on the edges. Their model is se- 
quentially updated, and for large values of the absorbing 
probability the system is rough (a « 0.5, and (3 « 0.33, 
compatible with their corresponding KPZ values), while 
for smaller values of the growing probability, the desorp- 
tion mechanism has a larger relative importance, small 
islands are more easily eliminated, and the asymptotic 
behavior is flat (a « /3 w 0). Moreover, these authors 
identify also a spontaneously broken symmetry in the 
model in the flat phase glj] . The phase transition in this 
case is also related to DP. 

More recently an apparently unrelated model has been 
formulated in order to describe fungal growth. A fungal 
colony grows invading an environment, with the peculiar- 
ity that the local growth probability is non-markovian 
^2| . With this main ingredient, and considering a trian- 
gular lattice, Lopez et al. |^^ found a roughening tran- 
sition between a rough EW phase and a flat regime. 

All the aforementioned models share a common prop- 
erty: The roughening transition can be related to a DP 
transition po|j2^ ]. Having introduced this class of DP 
related roughening transitions, in the next section we 
present a new model, in this same generic class, which 
turns out to be more suitable to be studied using the 
RSRG. 



III. THE MODEL 

Instead of trying to directly renormalize any of the 
previously defined models, we find it more convenient to 
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FIG. 1. System phase diagram. Filled circles (•) denote 
points in the Edward-Wilkinson universality class; crosses 
(x), points in the flat class; empty circles (o) points in the 
self-similar class and dashed lines points in the random depo- 
sition class. 

define a new model with the same generic phenomenol- 
ogy. The reason for this is that in the new model there 
is no proliferation of new parameters [Q when the RG 
transformation is applied, and therefore the RG flow and 
the fixed point structure can be analyzed in a simple 
way. This does not exclude, in principle, the possibility 
of renormalizing directly any of the previously described 
models. 

The model describes a surface driven by three dif- 
ferent physical processes: deposition, evaporation and 
condensation of particles on a one dimensional lattice, 
i = 1, . . . , L. At each lattice site we associate an inte- 
ger non-negative variable h(i). The random deposition 
process, corresponds to an external flux of particles, and 
occurs at each time step with probability p. The other 
two processes, i.e. evaporation and condensation, de- 
crease the height difference between neighbor sites, and 
constitute therefore a smoothening source. These occur 
with complementary probability 1 — p at each time step, 
and they are responsible for the generation of correlations 
among different sites. 

We consider flat initial conflgurations (hU) = Vi 
at i = 0), periodic boundary conditions [p4| (/i(l) = 
h{L) Vt), and sequential updating. The dynamics is de- 
fined by the following algorithm: At each simulation step 
a lattice site io is selected randomly. Its height can ei- 
ther be increased by one unit with probability p (random 
deposition) 



/i(io) -^ h{io) + 1 



(1) 



or alternatively, with complementary probability 1 — p 
the surface is smoothened away (evaporation or conden- 
sation) in the following way; 

/i(io) -^ h{io) + int[a (/i(io + 1) - /i(«o))] 

/i(io) ^ h{io) + int[a {h{io - 1) - /i(io))] (2) 

where a € [0, 1]. Each one of these two possible events 
occurs with equal probability (1 — p)/2. This process 
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FIG. 2. Evolution profiles for the different universality 
class at times (from bottom to top) t = 10, 10^, 10^, 10'' and 
10^ for L = 1000. (a) (p = 0.4, a = 0.5) corresponding to 
the Edwards- Wilkinson class (b) (p = 0.1, a = 0.5) in the flat 
class; (c) (p = 0.75, o = 1.0) in the self-similar class; and (d) 
(p = l,a = 0.5) random deposition class. Curves at t = 10"* 
and 10^ {t = 10, lO'^) are shifted upward (downward) for the 
sake of clarity. 



causes a decrease of the height difference (smoothening) 
between the site io and one of its neighbors. In the dy- 
namical rule given by Eq. (g) we have introduced the 
integer part function to enforce the height variables to 
take integer values. 

Let us now discuss the model phenomenology. We have 
investigated the parameter space (p, a) with analytical 
and computational methods in different points as shown 
in Fig.^ The evolution of the surface width can belong to 
one out of four different scaling regimes depending on the 
parameter values, p and a. These four universality classes 
define four qualitatively different growth morphologies, 
the corresponding typical surface profiles of which are 
shown in Fig. g. They correspond to the following uni- 
versality classes: The Edwards- Wilkinson |25|,|l|-p|, the 
flat, the self-similar, and the random deposition []l|-^ uni- 
versality class. 



A. Edwards- Wilkinson universality class 

The EW class is characterized by critical indexes a = 
1/2, P = 1/4 and z = 2 |]-|. We observe this type of 
scaling for parameter values p > Pc{a), (where p = Pc(o) 
is the separatrix between this phase and the flat one) 
with p ^ 1 and a ^ 0, 1 (see Fig. w where points in this 
class are marked with filled circles). In Fig. (a) we show 
the collapse of the surface width for different system sizes 
(see figure caption). The best curve collapse is obtained 
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FIG. 3. Width evolution for system sizes 

L = 32, 50, 64, 90, 128, and 256, for three different pairs of pa- 
rameter values corresponding to: (a) the Edwards- Wilkinson 
universality class, (p, a) — (0.4,0.5); (b) the flat universality 
class, {p,a) = (0.1,0.5); and (c) the self-similar universality 
class (p, a) = (0.75, 1.0). In (a) and (c) the width is rescaled 
with L" and plotted versus the time rescaled with L^, with 
a = 0.50 ± 0.01, /3 = 0.25 ± 0.01 in (a), and a = 1.00 ± 0.02, 
/3 = 0.50 ± 0.03 in (c). In (b) a semilogarithmic scale is used; 
there is a logarithmic dependence on time for t < 10, and 
a system size independent saturation width, compatible with 
a = 0.00 ± 0.01 and /3 ~ 0. 



for a = 0.50±0.01 and (3 = 0.25±0.01 in agreement with 
the system being in the Edwards- Wilkinson universahty 
class. This phase has a mean growing velocity different 
from zero (i.e. grows continuously in time) and has a 
rough surface appearance as shown in Fig. ||(a). 

Let us present here a simple argument showing why 
this rough phase is EW like and not KPZ like. For that, 
we calculate the mean over different runs of the quantity 
h{i,t + A) — h{i,t) (where Ai = 1/i is a time increment). 
We can write: 

{h{i,t + At) - h{i,t)) = 
P , 1 (1-p) 



L L 2 
1 (1-p) 



a {h{i + l,t) - h{i,t)) + 



L 2 
L 2L 



a {h{i- l,t)-h{i,t)) 



(3) 



(4) 



a {h{i + 1, i) + h{i - 1, i) - 2h{i, t)) (5) 



where it can be seen that there are two different contribu- 
tions to the local velocity at each site i; the contribution 
coming from deposition, proportional to p, and the term 
arising from the smooth away processes. Observe that in 
this calculation we have neglected the effect of the integer 
part function, which we assume to be irrelevant in this 
phase (i.e. we expect it to reduce the mean velocity but 
not to change qualitatively the behavior). 

Dividing Eq. (0) by At and performing the thermody- 
namic limit, i.e., L -^ oo (or equivalently At -^ 0), we 
obtain: 



dh 



(1-p) 



a V^/i. 



(6) 



that is the Edwards- Wilkinson equation driven by an ex- 
ternal force p. For p = 1 or a = this equation describes 
random deposition as, in fact, is the case in our model. 
Eq. (0) fails to describe the growth of our model in the 
cases p < Pc OT a = 1. In fact, as stated before, in Eq. 
(^ the integer part function has been neglected and this 
approximation is incorrect when: (i) p < Pc, the fluc- 
tuations of the surface are small, and the rounding off 
mechanism due to the integer part function takes over, 
pinning and flattening the surface, and (ii) in the case 
a = 1, for which, in the absence of smoothening, the 
continuum limit has to be taken more carefully. In this 
last case, we expect the dynamics to be controlled by the 
diffusion of height steps. 



B. Flat universality class 

This universality class is characterized by critical in- 
dexes a = and /? = (with logarithmic corrections). 
We observe this phase for parameter values: p < Pc(a) 



with p ^ and a ^ 0, 1. In Fig. |] we plot points in this 
class with crosses in the parameter space. In this phase 
height fluctuations are independent on sample size (i.e. 
a = 0), the mean surface velocity is zero in the thermo- 
dynamic limit, and the lowest level {h — 0) is occupied 
with a finite density as shown in Fig. 0(b). Due to fi- 
nite size effects, finite systems in this class may have a 
non- vanishing velocity. 

In Fig. 0(b) we have plotted in a scmilogarithmic scale 
the surface width versus time for different system sizes 
(see figure caption). Observe that the width has a loga- 
rithmic dependence on t for short times 



W{t) ex log(t). 



(7) 



For all the different points in this phase, by employing 
data collapse techniques we evaluate a = 0.00 ± 0.05 and 
/3 = with logarithmic corrections. 

The existence of this phase is due to the fact that for 
small values of p the deposition process is much unlikely 
than the smoothening one, and there is a physical con- 
straint preventing the surface to go below the lowest, 
h — Q level. The integer part function in the dynamic 
rules favors the evaporation process by giving and extra 
negative drift term with respect to Eq. ||. This extra 
term binds the surface to the lowest level making it flat 
for small values of p. In particular, for a fixed value of a, 
the roughening transition corresponds to the value of p 
for which this extra negative term is equal to —p. For val- 
ues of p slightly larger than that critical value the surface 
unbinds and grows with constant velocity. 



C. Self-similar universality class 

This class is characterized by a = 1, (3 — 1/2 and z = 
2, it is observed for parameter values a — \ and p ^ 0, 1 
which are plotted in Fig. [| with empty circles. In Fig. 
H(a) we show the collapse of the width for different system 
sizes (see figure caption); from the best data collapse we 
measure a = 1.00 ± 0.02 and (3 = 0.50 ± 0.03. The 
surface in this regime is therefore self-similar, that is, 
height fiuctuations are of the order of magnitude of the 
system size (see also Fig. ||(c) where a typical profile is 
shown). In this case, contrarily to a ^ 1 we observe no 
roughening transition (except for a trivial one at p = 0.) 

Observe that for a = 1 the dynamical rules can be 
written as: 

h{if) , t) = h(io , t — 1) + 1 with prob. p 

/i(io ,i) = /i(io + l ,i- 1) with prob. (l-p)/2 (8) 

h{io ,t) = h{io-l ,t-l) with prob. (l-p)/2, 

where the integer part function, being redundant, has 
been omitted. Therefore, the extra negative driving 



term, arising as a consequence of the integer part func- 
tion in the dynamics (which, as discussed in the previous 
subsection favors evaporation) does not exist in this case. 
This explains why there is no roughening transition for 
a = 1. 

The main quaHtative physical difference between this 
phase and the EW and the flat phases is that here there 
is no smoothening of height gradients. Large steps diffuse 
in space due to the effect of the second and third rules 
in Eq. (ph, but do not smoothcn away, creating a much 
rougher surface. 

Now we study the mean field solution of this case, ob- 
tained by neglecting spatial correlations. This turns out 
to be an useful calculation as we will show afterwards. 

The master equation for the probability P{h,t) for a 
given site to have a height value h at time t is 



dP{h,t) 
dt 



= pP{h-l,t)-pP{h,t)+ 

+ {l-p)P{h,t)~{l~p)P{h,t) 
= p[P{h-l,t)-Pih,t)] 



(9) 



where we have assumed nearest neighbor sites to have 
the same probability distribution as the site under con- 
sideration. Eq. (|9|) is the same equation that we would 
obtain if we considered a mean field approximation of 
the random deposition process. In order to solve (g), it 
is convenient to introduce the generating function defined 
as 



G{x,t) = ^x^P{h,t) 



h=0 



with the normalization condition G(l,i) = 1 Vi. 
master equation written in terms of G{x, t) reads 



(10) 



The 



dG{x) 
dt 



= p [xG{x, t) - G(x, t)] = p{x - l)G{x, t). (11) 



By taking derivatives with respect to the dummy variable 
x, one can write 



W^it) = <h^-<h>^ > = 

= [dlGix)+d,Gix)]\^^^-[d,Gix)l__ 



(12) 



Taking time derivatives of both sides and integrating the 
resulting equation, we finally obtain W(t) = \ppi^ and 
consequently /3 — 1/2. On the other hand, the height 
differences among neighbors (or steps) perform a random 
diffusion in the direction perpendicular to the growth. 
Like in the case of a random walker the average step dis- 
placement scales as square root of time. But the effective 
time for step displacement is (1 — p)i, and therefore in 
order to cover a distance L, a characteristic time 
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FIG. 4. Saturated width of a system of size L = 64 with 
a — 0.5, as a function of p. The data are extremely well fitted 
by the mean field prediction 



is needed. Consequently the dynamical exponent is z = 
2. Combining this result with the fact that p = 1/2 we 
obtain a = 1 just by using the scaling relation a — z/ j3. 
Therefore, the values of the critical exponents found in 
mean field approximation agree perfectly with the nu- 
merical results. Moreover, the mean field approximation 
also reproduces the dependence of the saturation width 
as a function of p for a fixed system size. In fact, using 
W(t) = ^/pt and Eq. (|l|) we find: 



Wsat{p) 



P 



(1-p) 



(14) 



(13) 



in perfect agreement width the simulation results as 
shown in Fig. 0. 



D. Random deposition universality class 

The random deposition class is characterized by a = oo 
and (3 = 1/2. With dashed lines in Fig. |l|we have plotted 
the two lines in this class, namely those given by the con- 
ditions a = and p= 1 respectively. For these values the 
only physical mechanism present is random deposition of 
particles. In this phase the surface is spatially uncorre- 
lated, and its roughness increases rapidly with time, with 
f3 — 1/2, and there is no width saturation (z = oo) as 
shown in Fig. ||(d). 



E. The roughening transition 

As stated previously, both in our model and in the pre- 
viously discussed ones the transition between the flat and 
the rough phase belongs to the directed percolation (DP) 
universality class. The reason for this can be explained 
by mapping the dynamics in a directed percolation like 



model. In order to do this in our model (for the others 
similar mappings can be performed), we introduce a new 
set of variables {si} defined as 



s, == 1 if h{i) = 
Sj = if h{i) ^ 0. 



(15) 



By studying the evolution of these variables we focus our 
attention on the lowest height level occupation. This in- 
formation is enough to describe the phase transition p8[ , 
as we show in what follows. A site i with Si = 1 corre- 
sponds to an occupied (active) site in a DP like model (or 
better, in a contact process model which is a sequentially 
updated version of DP f^y]). On the contrary, Si — for 
empty (absorbing) sites of the DP like model. The depo- 
sition process may change with probability p an occupied. 
Si = 1, site into an unoccupied, Si = 0, site. This same 
mechanism is also present in the contact process [|20|] . On 
the other hand, the smoothcning mechanism, occurring 
with probability 1 — p, may induce Ig^l the " infection" 
of an empty site by a neighboring active site as in the 
contact process. In this language, the flat surface phase 
corresponds to the DP active phase, and the rough phase 
to the absorbing one. The key feature is that in absorbing 
regions (this is, in regions with Si = 0), activity cannot 
be generated spontaneously. This is the main physics of 
DP and therefore the critical behavior at the roughening 
transition is related to the transition into an absorbing 
phase of DP and related models. 

By using this mapping the scaling of some magnitudes 
can be related to DP exponents. In particular, the den- 
sity of sites at the lowest level, no, should scale as a 
function of the distance to the critical point, e = \p — pc\i 
like no ~ e^°^ . Right at the critical point no decays in 
time as no(t) ^ i~^°^. For finite systems, of size L, 



no ^ L 



-PdP /l'±,DP 



L- 



(16) 



where i'±^dp is the correlation length exponent. Finally, 
the mean surface velocity is inversely proportional to the 



life time r 



II, -DP 



of the DP active phase, where 



^11, DP is the usual correlation time exponent. Therefore, 

'"11 ^ T^-l'll/'^-L 



V ' 



L- 



(17) 



The one-dimensional DP values of the previously intro- 
duced exponents are: (3dp = 0.27649(4), = 0.15947(3), 
Xf = 0.25208(5) and x^ = 1.58074(4) (numbers in paren- 
theses denote uncertainties in the last figure) [ pTj . 

In order to verify numerically if the previous predic- 
tion holds we have studied extensively the a = 0.5 case. 
We have performed numerical simulations in order to de- 
termine critical exponents, to be compared with their 
corresponding DP values. In particular, by considering 
a system size L — 3000, and averaging over 8000 in- 
dependent runs, we measure no as a function of time. 
A power law decay is observed at (the critical point) 
Pc = 0.13740 ± 0.00005, with an associated exponent 
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FIG. 5. Critical indexes characterizing the DP transition. 
(a) Log-log plot of the density of sites at lowest level as a 
function of time at (0.1374, 0.5); the associated critical expo- 
nent is 6' = 0.160 ± 0.003. (b) Log-log plot of the asymptotic 
density of states at the lowest site versus e = Pc — p\ from 
the slope we measure Pdp = 0.275 ± 0.007. In (c) and (d) we 
report the density of sites at the lowest level and the mean 
velocity, at the critical point for different system sizes. The 
corresponding critical indexes are Xf — 0.250 ± 0.003 and 
Xv = 1.57 ±0.01. 



= 0.160 ±0.003 



(18) 



in perfect agreement with its DP value. Using the pre- 
viously estimated value of the critical point, we mea- 
sure the stationary density of sites at the lowest layer 
h = ior p < Pc- In Fig. H(b) we show the data for 
e € [10~*, 10~^]. Performing a fit for small values of e 
we find 



^DP = 0.275 ±0.007 



(19) 



also in good agreement with its DP value. 

By performing a finite size scaling analysis, we de- 
termine Xf = 0.250 ± 0.03 and x„ = 1.57 ± 0.01 to 
be compared with their corresponding DP values Xf = 
0.25208(3) and Xy = 1.5807(5) respectively. Numerical 
data are showed in Fig. S (c)-(d). 



Summing up, all the measured exponents confirm the 
hypothesis that the roughening transition is controlled 
by a DP fixed point. 



IV. RENORMALIZATION APPROACH 

Having described in detail the model phenomenology, 
in this section we apply the RSRG technique recently 
proposed |12-14| to renormalize our model. For a de- 
tailed description of the method we refer the reader to 
reference 14 1. In a nutshell the RG method consist of 
two main ingredients. One is the definition of cells of 
generic substrate length L^ = 2^Lq (where Lq is the 



length of the minimal relevant substrate scale) and height 
/ife, with which the space in which the surface grows can 
be covered; cells are progressively invaded by the grow- 
ing surface. The second one is the identification of the 
effective dynamics of these cells at a generic scale; this is, 
one studies the effective rules by which cells are progres- 
sively invaded by the growing surface at a given scale. 
This effective dynamics is described by parameters that 
change upon changing the scale of description. In our 
case we have a set of parameters x^ = (afe,Pfe) at scale 
Lfc. Known the effective dynamics at scale Lfe, defined by 
Xfe and the associated width at that scale Wk , the renor- 
malization problem consists in determining what are the 
width values Wk+m and the effective dynamic parame- 
ters Tik+m at a coarser scale, Lk+m- This problem was 
analyzed in |13-n3, with the final conclusion that one 
can write 



W^ 



W?F„ 



(Xfe) 



(20) 



where Fm(xfe) = 1-1- 4M;^(m,Xfc) and w{m.,Xk) is the 
width of a system composed oi 2"^ cells of unitary height 
with a dynamics driven by the parameters x^ . This last 
functions can be determined by means of a rather inex- 
pensive MonteCarlo simulation. The width at scale k + 2 
can be calculated in two different ways; (i) by direct use 
of the previous formula with m = 2, or alternatively, 
(ii) by iterating twice the previous transformation with 
m = I. Imposing that both the previous procedures give 
the same result for the width Wk+2 , one obtains a renor- 
malization condition, namely, f|l2[|l4 



Xk+l 



= F^'[F2{xk)/Fi{xk)]. 



(21) 



For mono-parametric models this equation is enough to 
determine the evolution under RG transformations of the 
effective parameter Xk , and from it all the scale invariant 
physics can be elucidated |12|-p^. In the present two- 
parametric dynamics it is necessary to consider another 
independent analogous equation, corresponding to calcu- 
lating the width at scale A; -|- 3 in two different ways; 



Xfc+l 
Xfe+l 



rM^2(xfc)/Fi(xfc)] 

2-i[F3(xfc)/Fi(xfc)]. 




FIG. 6. A system of size Lk+2 = 4I/fe can be studied as 
composed of four cell of size Lt , or alternatively by iterating 
twice a partition in cells of half length. See text. 

The first equation is just Eq. (pT|), and the second states 
that the width obtained by dividing the system in 8 cells, 
should be the same as the width obtained dividing first 
the system in two blocks, and then each of these blocks 
in four sub-cells. This set of equations give the (discrete) 
flux of the effective dynamics parameters upon coarse 
graining. The scale invariant dynamics is determined by 
the fixed point (or points) of Eq. (p2|), x* for which 



Fi(x*) = F2(x*)/i^i(x*) 
i^2(x*) = F3(x*)/^i(x*) 



(23) 



is satisfied. Once x* is known the a exponent is easily 
determined by 



«-^log2(^l(x*)). 



(24) 



(22) 



In order to evaluate the Fi, F2 and F3 functions we 
have performed MonteCarlo simulations for system of 
size L — 2,4,8, and averaged, once the stationary state 
is reached, for long enough times. We have determined 
these functions on the sites of a 100 x 100 lattice in the 
parameter space, this is, for values of p and a multiples 
of 0.01. 

In Fig. we show the curves corresponding to Fi — 
F2/F1 and F2 — F3/F1. The intersection points of these 



0.0 0.2 




r,=F7r, 

f. (0.75,1.00) 



(0.14,0.50) (0.40,0.50) 



FIG. 7. The two lines (marked respectively with dia- 
monds and empty circles) correspond to the values for which 
_Fi(x) = i^2(x)/7^i(x) and F2(x) = f3(x)/Fi(x), within the 
numerical error respectively. The intersection points of these 
two lines define the dynamic fixed points, invariant under RG 
fiow. See text. 

curves are the RG fixed points within numerical accuracy. 
We find the following stable fixed points: 



p = 0.75, a = 
p = 0.40, a = 
p = 0, Va 



1 corresponding to a 

0.5 corresponding to a 

corresponding to a 



1.00 ±0.02 
0.50 ±0.01 (25) 
0.00 ±0.01 



while at (p, a) = (0.14,0.5) there is an unstable fixed 
point. The stability is determined by observing whether 
nearby points fiow into the fixed point or fiow away from 
it, under application of Eq. (E2h, and the corresponding 
values of a are determined using Eq. (p4|). 

We have also determined different rcnormalization flux 
lines as shown in Fig. @. Observe that each continu- 
ous flow line in Fig. g has been constructed by joining 
together consecutive points obtained from the rcnormal- 
ization group equations |2j (which determines a discrete 
iterative mapping and not a continuous fiow). Observe 
that the lines defined by a = 0.5 and a = 1 define invari- 
ant manifolds. 

The accordance between the rcnormalization fiow di- 
agram (Fig. 0) and the phase diagram independently 
found with simulation on large systems (Fig. M is very 
good (compare Fig. ^and^). In particular: 

• For all points yielding in the EW phase we find 
trajectories converging to the stable fixed point 
(0.4, 0.5) with a = 0.5 ± 0.01, in perfect agreement 
with the EW value. 

• The separatrix in Fig. |l| is reproduced quite accu- 
rately in Fig. y, in particular, the unstable fixed 
point (0.14,0.5) is located on it. 

• Points to the left of the separatrix (i.e. points in 
the fiat phase) fiow towards the fixed line at p = 
for which a = (as corresponds to the fiat phase). 




FIG. 8. Renormalization Group fiow for Eq. (0). The 
stable fixed points are p — 0.4, a = 0.5 (EW behavior); 
p = 0.75, a — 1 (self-similar phase); p,Va (flat phase); and 
the unstable point p — 0.14, a — 0.5 (corresponding to a DP 
transition). For the fiux line starting at (0.95, 0.95), we show 
explicitly the points found by iterating the renormalization 
group equations. To guide the eye discrete sequences of points 
obtained by iteration of the RG transformation are plotted as 
continuous lines. 



• All points with a = 1 and p > 0, i.e. in the self- 
similar phase, fiow to a fixed point at (0.75, 1) with 
a = 1.00 ±0.01. 

We have verified all the above conclusions to be sta- 
ble when more refined RSRG algorithms arc considered. 
This is, instead of considering by-partitions and quadri- 
partitions of a given growing surface, one can consider 
larger partitions, and the technique described above re- 
mains the same in spirit [n2Ul4|. By doing this, we ob- 
serve all the above described results to remain unaltered. 
This stability upon changes in the RSRG details supports 
the fact that no new parameter has to be introduced upon 
renormalization to describe he scale invariant dynamics 
p4| , i.e. the scale invariant dynamics is well described at 
an arbitrary scale by the found fixed points. 



V. CONCLUSIONS 

We have discussed different growth models exhibiting a 
one-dimensional roughening transition. In particular, we 
have introduced a new two-parametric model capturing 
the main physics of this type of roughening transition, 
i.e exhibiting a rough and a flat phase, and a critical line 
separating them, related to directed percolation. The 
new model has the technical advantage that no new pa- 
rameter proliferate upon renormalization. Apart from 
the previously discussed phenomenology, the new model 



exhibits also a self-similar phase. All the different phases 
and transitions have been analyzed by means of extensive 
numerical simulations and some analytical approaches. 

We have renormalized this model by using a recently 
introduced RSRG approach. In particular, fixed points, 
corresponding to the scale invariant dynamics are found, 
and the corresponding roughness exponent determined. 
The results are in perfect accordance with all the nu- 
merical and analytical findings. In particular, the phase 
diagram is perfectly reproduced: The rough, flat, and 
self-similar phases, as well as the separatrix among them, 
and their associated roughness exponents are identified 
with great accuracy. 

This confirms the general validity of the RSRG method 
to deal with anisotropic fractal growth in cases others 
than KPZ growth for which it was explicitly design. 
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contract FMRXCT980183. 
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